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Abstract 

This paper describes a quantum algorithm for finding the maximum 
among N items. The classical method for the same problem takes 
O(N) steps because we need to compare two numbers in one step. This 
algorithm takes 0(>/~N) steps by exploiting the property of quantum 
states to exist in a superposition of states and hence performing an 
operation on a number of elements in one go. A tight upper bound of 
for the number of steps needed using this algorithm was found. 
These steps are the number of queries made to the oracle. 



1 Introduction 

Quantum computing as a new paradigm of computing was introduced in the 
early 1980's in the work of Feynman Q and Beinoff pL Later they were 
shown to be atleast as powerful as classical computers. But the real inter- 
est in exploring the bridge between physics and computation arose when 
quantum algorithms which improved over their classical counterparts were 
proposed. One was the algorithm for factorization in polynomial number 
of steps by Shor || and another was algorithm which achieved quadratic 
speedup for the classic problem of database search by Grover 0|. This 
power of quantum computation is due to the fact that the state of a quan- 
tum computer can be a superposition of basis states and we can perform 
a operation on multiple quantum states simultaneously. This paper applies 
the quantum techniques to a very basic problem in order statistics i.e. of 
finding the extrema of a set of numbers. 
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2 Problem 



Let T[0 . . . N — 1] be an unsorted table of N distinct items. We have to find 
the index i of the element such that T[i] is maximum. 

3 Algorithm 

The algorithm below uses an oracle that computes: 




1 if T[j] > T[i] 
otherwise 



This algorithm starts out with an initial guess of the index of the maximum 
and using Grover's search algorithm, finds out a new element which is the 
index of one of the elements greater than the initial indexed element. The 
new index forms our new guess and the whole procedure is repeated. The 
various steps of the algorithm are : 

1. Choose a index y randomly from the set {0..iV — 1}. This forms our 
initial guess of the index of the maximum element. 

2. Repeat the following 0(\^N) times: 

• Initialize the state as \if>) = J2i 77^ I*) \u)- This can be achieved 
in 0(logN) steps using a n-bit Hadamard transformation. 

• Apply the Grover's search algorithm for finding a marked state 
(if no. of marked states is > 1) We consider a state x as 
marked if the following is satisfied, 

fy(x) = 1- 

• Make a measurement on first register. Let the result of the mea- 
surement be xq. Then our new guess y will be xo- 

3. Return y as the index of the maximum element in the table. 

4 Complexity Analysis 

Let E(N, t) be the expectation value of the number of iterations for finding 
the maximum of ./V items in which t are marked. Then we write a recurrence 
equation for E(N, t) as : 
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1 N 
E(N,t) = - [E(N, t — 1) + ....E(N, 1)] +6W- (1) 

(Since the expected no. of iterations to find the index of a marked item 
among N items where t are marked is at most 6^/N/t) |5[. Writing the 
equations when the no. of marked items is t and (t — 1), we get 

E(N,t) = - t Y,E(N,*)+^j (2) 

i=l ' 



1 s / N 



E(N, t - 1) = — Y, E{N, i) + 6W — (3) 



Multiplying (2) by t and (3) by (t — 1) we get, 



t-i 



tE(N,t) (4) 



t-2 



(t - 1)£(JV, 4-1) =S ^ + 6 V^(* " 1) (5) 



i=l 



Subtracting (5) from (4) and rearranging, we get 



£7(JV, t) =E(N,t-l) + 6^-(Vi- Vt^l) (6) 

Writing the same equation for (t — 1 ),..., 2 and adding all of them, we get, 

E(N,t) = E{N,l)+6^Ny-(Vi-Vi^l) (7) 
This can be rewritten as : 



E(N,t) = E(N,l)+6VNj2-^(l- \A--) (8) 



i=2 



Applying the relation (1 — x) n > (1 — nx) and replacing summation by 
integral, we get 

E(N,t) < E(N,1) +6\/iV / =cfe. (9) 

Ji=i 2ztJz 



3 



Evaluating the integral, we get 

E(N,t)<E(N,l) + 6y/N(l-^p). (10) 

yt 

An upper bound is achieved for t = N — 1, when we get 

E(N,t) < 6.8v^N 

Thus, we get the upper bound on the number of iterations as 6.8\/iV. After 
we had the initial idea of this algorithm, Dr.Grover told us that there already 
exists a minimum searching algorithm by Diirr and H0yer || . We found that 
both the algorithms were the same. Our complexity analysis yields a better 
bound than the original one of 15yN. We can repeat this algorithm to 
obtain the maximum with a probability of 1 — (^) fc after k repetitions of this 
algorithm. 

We next show a high probability bound for the problem. Applying 
Markov's inequality, we have, 

P{X > k{X)) < ~ 

where {X) is the expected value of the random variable X. Therefore using 
the upper bound on the expected number of iterations, we get 

P{X > 6.8kVN) < - 
k 

So after 13. steps, we would have reduced the probability of failure to 
less than 1/2. Repeated iterations of this algorithm give us a high probability 
bound but with an additional factor of 0(logN). It would be of interest to 
improve this bound further. 
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